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We point out that stars in the mass window ~ 8 — 12Mq can serve as sensitive probes of the axion- 
photon interaction, QA-yy- Specifically, for these stars axion energy losses from the helium- burning 
core would shorten and eventually eliminate the blue loop phase of the evolution. This would 
contradict observational data, since the blue loops are required, e.g., to account for the existence 
of Cepheid stars. Using the MESA stellar evolution code, modified to include the extra cooling, we 
conservatively find <5U 77 < 0.8 x 10~ 10 GeV -1 , which compares favorably with the existing bounds. 
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INTRODUCTION 

For a particle physicist, stellar interiors represent ex- 
tremely hermetic detectors, sensitive to very rare pro- 
cesses. For example, the 7* — > vv decay [1] measur- 
ably drains energy from the core of a red giant star, even 
though the probability of this decay to occur between 
successive interactions of the plasmon 7* is only ~ 10 -26 . 
Importantly, the rate of the energy drain cannot be too 
different from the Standard Model (SM) prediction, al- 
lowing one to constrain neutrino properties [SHI] (see [5] 
for further references). The same argument extends to 
new physics scenarios with light, weakly interacting par- 
ticles [5J. Numerous examples include majorons [7], light 
supersymmetric partners [SJ, novel baryonic or leptonic 
forces [9], and more recently unparticles [10] and extra- 
dimensional photons [11] . A particularly compelling sce- 
nario is furnished by the axion |121 113) , which is the sub- 
ject of this letter. Below, we reexamine the astrophysical 
implications of the axion-photon coupling and point out 
that, contrary to the standard lore, stars with masses 
r-j 8 — 12M Q are very sensitive to it. 

The axion arose from a proposal to account for the ab- 
sence of CP violation in the strong interactions (QCD) 
[T3HTT] . The SM QCD Lagrangian admits a CP-violating 
GG term, which, if present, would impact physical ampli- 
tudes through nonperturbative effects [18T[2Tj . In partic- 
ular, one may expect the neutron to have a large electric 
dipole moment [3H [53], contrary to observations [23] , 
The axion proposal addresses this by promoting the co- 
efficient of the GG term to a dynamic field, which is con- 
structed to be the Goldstone component of a U(l) field. 
The corresponding broken symmetry (Peccei-Quinn) is 
anomalous, hence the Goldstone couples to the SM fields, 
particularly the pion, and gains a small potential. This 
potential dynamically drives the axion field to the CP 
conserving vacuum, solving the problem. 

Being a pseudo- Goldstone boson, the axion can be light 
enough to be produced in stars. More precisely, the axion 
mass vtiA and decay constant Ja are related to those of 
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FIG. 1: The Figure shows why He burning is optimal for 
probing the A"/"/ coupling. The curves trace the evolution of 
the central density pc and temperature Tc in stars of 1 Mq 
(bottom, red) and 15 Mq (top, blue). Different burning stages 
are labeled. The shaded region shows the range of conditions, 
for which axion emission with gio = 1 contributes at least 
90% of the non-photon energy loss. The dashed isocontours 
correspond to loss rates of 10 1 and 10 4 erg/g/s, as labeled. 

the pion, as to^/a — m„f n , or [TBI |2"5I |2"6] 

(m A /leV)(/ j4 /10 7 GeV)~0.6. (1) 

Below, we will be interested in axion emission from He- 
burning stellar cores, which have temperatures ~ 10 4 eV. 
Eq. ([T]) then tells us that for Ja above the weak scale the 
axion is indeed light enough to be thermally produced. 

The high scale of Ja also ensures the second condi- 
tion: axions interact weakly enough to free-stream out 
of stellar cores. The couplings of the axion field 4>A 
to axial SM currents = ^> fj^j^^ /, and to pho- 
tons are both suppressed by Ja- £ € C//J Jfd^cpA + 
C 1 a/(8irfA)4>AF fJ ,vF f2 ". It is easy to verify that axions 
emitted from _ffe-burning stellar cores do not reinteract. 

In this letter, we specialize to the second coupling, 
A77. In many motivated axion models the dimension- 
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FIG. 2: The evolution on the Hertzsprung-Russell (HR) diagram of stars with masses 7 — 12 Mo, with standard cooling (left) 
and with the addition of the ^77 coupling at the edge of the sensitivity of CAST (right). For stars 9Mq < M < 12M Q the 
axion losses completely eliminate the blue loop stage. The dashed lines show the instability strip (conditions for Cepheids). 



less coefficient C 7 is 0(1). For example, for the well- 
known KSVZ [23 US] and DFSZ [30] scenarios, we 
have \C y \ ~ 1.9 and 0.7 respectively. It follows that the 
mass of the axion and the coupling strength to photons 
Ga-ty = Ct-cV^tt/a) are proportional to each other. For 
example, for the KSVZ model one has 

G^ 2 ~ 3.7 x 1(T 10 GeV(m A /l cV). (2) 

For large Ja the axion becomes very difficult to de- 
tect ( "invisible" ) and the coupling Ga 77 becomes one of 
the key experimental handles [31] . Many recent searches 
have been targeting this coupling 32 [59] , including dark 
matter detectors, DAMA [33] and CDMS [3J], dedicated 
axion telescopes, Tokyo [35] and CAST 36J, a reactor 
experiment, TEXONO [37] . and even a solar neutrino 
experiment, Borexino |38| . These experiments exclude 
various segments on the line in Eq. ^ in the range 
10° < itia < 10 5 eV, corresponding to the Peccei-Quinn 
scale 10 2 GeV < f A < 10 7 GeV. 

Remarkably, a stellar cooling bound published over 20 
years ago |39j excludes this entire range, pushing the 
bound on Ja all the way up to 10 7 GeV. Given the ob- 
vious importance of this result, it is highly desirable to 
confirm it with more than one type of stellar systems. 
This is done below. The bound of [31] makes use of low 
mass (~ I.3M0) stars. We show that stars of heavier 
masses, 8 — 12M Q , can also be used as effective axion 
probes, an rather unexpected result (c/. [40], p. 37). In 
these stars, axion cooling can qualitatively change the 
evolution, with clear observational consequences. 

We strengthen the astrophysical bound also in another 
way. For our numerical modeling, we make use of the 
publicly available and community tested MESA stellar 
evolution code, to which we release our modifications cap- 
turing the axion cooling rates. Our analysis can thus be 
independently verified and - we hope - extended. 



WHY HELIUM BURNING? 

The axion-photon coupling leads to energy loss via the 
Primakoff conversion [T3 HT] : photons convert into ax- 
ions in the background field of nuclei. The conversion 
rate is controlled by the finite range of the Coulomb field 
in plasma, which regulates what would otherwise be a 
forward scattering logarithmic divergence |42| . The re- 
sulting expression is well established [S] 02] ; in a nonde- 
generate medium, per unit mass, the axion loss is 

e A = Z(e)%?^ = 27.2^- s Z(e)gl Tlp- 3 \ (3) 

where g w = G A77 /(10 10 GeV" 1 ), p 3 = p/(10 3 g/cm 3 ), 
Tg, = T/10 8 K. Three powers of temperature come from 
the photon number density, one from the energy loss per 
photon, and the remainder from the form of the (plasma- 
regulated) cross section. 

The coefficient Z(£ 2 ) is a function of £ 2 = (k s /2T) 2 , 
with Ks being the Debye-Huckel screening wavenumber. 
Z(£ 2 ) is given as an integral over the photon distribution 
(see Eq. (4.79) in [S]) and is generally C(l) for relevant 
stellar conditions. For example: for the Sun, £ 2 ~ 12 
and Z ~ 6; for the low-mass He burning stars, £ 2 ~ 2.5 
and Z ~ 3 [5]; finally, for a 1OM He burning star of 
interest here, £ 2 ~ 0.1 and Z ~ 0.4. To include the axion 
losses in the stellar evolution code, we need a simple, 
yet accurate, parameterization for Z(£ 2 ). Obviously, this 
function needs to interpolate between the limits Z(£ — > 
0) = (7r 3 /30)£ 2 ln(3.99/£ 2 ) and Z(£ 2 -> 00) = 2tt 5 /63, 
but the interpolation needs to also accurately reproduce 
the intermediate regime, since the physically interesting 
values of £ 2 lie there. We propose using 

ztt 2 ) - (nrrife + aroffe) - + W) ■ ( 4 ) 
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The accuracy of this parameterization is better than 2% 
over the entire range of £. 

Using the cooling rate in Eqs. (3|4|, we plot in Fig. [I] 
the region where the axion cooling with gio = 1 com- 
prises at least 90% of the overall non-photon energy loss. 
The effect of the axion is pronounced at moderate tem- 
peratures and densities, ordinarily the domain of photo- 
production (je" —¥ e~vv)\ for higher temperatures, it is 
overtaken by the SM pair production (e + e _ — > vv), while 
for higher densities, the SM plasmon decay dominates (c/. 
[321 H3] ) . Since the rate of axion emission increases with 
temperature (as illustrated by the two dashed isocon- 
tours), the optimal temperatures for probing Primakoff 
losses are generally in the upper part of the shaded re- 
gion, 1 x 10 8 K < T c < 4 x 10 8 K. These are precisely 
the conditions at which Helium burns. For illustration, 
we show two curves depicting the evolution of the cen- 
tral temperature and density in 1 M Q and 15 M Q stars. 
The calculations were carried out with the MESA code, 
without the axion cooling. Next, we show what happens 
as this cooling is added. 



81). Notice that in this case there is no qualitative change 
of the evolution. The argument constraining the axion is 
based on counts of low-mass HB stars in stellar clusters 
and in the galactic disk. Faster burning of He due to 
axion losses would reduce the counts. 

We next turn to our main calculation, the impact of 
the axion on stars of masses ~ 7 — 12Mq. The results are 
shown in Fig. [2] as the evolutionary tracks for these stars 
in the Hertzsprung- Russell (HR) diagram for gio = and 
9io = 0.88. The second value represents the limit of sensi- 
tivity of the CAST experiment for very small axion mass 
(off the KSVZ line). We see that even such small axion 
coupling qualitatively changes the evolution. Normally, 
these stars, after reaching the red giant tip, travel back 
to the left (blue) side of the HR diagram. This is the 
well-known blue loop phenomenon 50 54 . With the ax- 
ion cooling, however, for stars with 9M Q < M < 12M@, 
this evolutionary stage disappears altogether. 



DISCUSSION 



AXION COOLING IN MESA EVOLUTION CODE 

MESA ("Modules for Experiments in Stellar Astro- 
physics") is a robust, open source, modular ID stellar 
evolution package [33] • Its release several years ago rep- 
resents a very significant development in the field of stel- 
lar astrophysics. The "instrument paper" [45] has shown 
MESA to be remarkably versatile, capable of handling 
not only solar-mass stars, but also objects as diverse as 
10- 2 Af o and 1O 3 M . 

MESA has been widely accepted by the astrophysics 
community. Yet, its applications to studying new physics 
have so far been confined to models of modified grav- 
ity [35 -48 . To the best of our knowledge, the present 
letter represents the first use of MESA to constraining 
new particles. Given its openness and robustness, we are 
convinced that over time MESA will become a standard 
tool for probing fundamental physical processes in stars. 

We added the axion cooling effect, Eqs. ( 3p l to the 
MESA code (specifically, to the neutrino energy loss rou- 
tine neu.f). Our modified code is being made publicly 
available [3S], so that our analysis can be reproduced 
and further extended. As a simple verification test, we 
ran a 1.3 M Q model with and without the axion cool- 
ing, paying particular attention to the duration of the 
He burning (the Horizontal Branch, HB, stage). This is 
the model studied in [35] and which has been the basis 
for the published bounds over the last two decades. The 
model ran straightforwardly through Hydrogen burning, 
Helium flash, and the HB stage. The resulting durations 
of the HB stage were in excellent agreement with [3"9"] : 
1.2 x 10 s yrs without the axion 0.7 x 10 8 yrs upon adding 
axion-photon coupling with gio = 1 (see also [3H], page 



Let us examine the physics behind the disappearance 
of the blue loop. It is helpful to look at the evolution 
of the internal structure through the He burning stage. 
Fig. [3j shows this evolution for a representative 9.5M Q 
star, with the horizontal axis showing the stellar age in 
millions of years (Myr). The top panel corresponds to 
the standard case (gio = 0), while the middle and the 
bottom one have gio = 0.6 and 0.8 respectively. 

In all three models, Helium is ignited at 23.8 Myr and 
the preceding evolution is not noticeably changed by ax- 
ion losses. The axion losses, however, do speed up the 
He burning stage, as expected (more losses require faster 
burning). The duration of this stage decreases from ~ 2.6 
Myr for 510 = to - 2 Myr for g w = 0.8. Notice that 
in the top panel, at 25.5 Myr the star undergoes a trans- 
formation: its radius contracts, while its surface tem- 
perature rises. This is the blue loop phenomenon: the 
star transitions from a red giant with a large convective 
envelop to a more compact blue giant with a radiative 
envelop. The same transition also occurs in the mid- 
dle panel, but notice that the modest overall shortening 
of the He burning stage significantly shrinks the blue 
loop stage. In the bottom panel, the He burning stage 
is shortened enough that the blue loop simply does not 
have time to start. The core exhausts its Helium, then 
contracts to Carbon ignition (cf. Fig. [IJ, at which point 
SM neutrino losses increase so much that the rest of the 
evolution proceeds in a very short time (see, e.g., [55 ). 

This is the basis of our argument for constraining the 
axion-photon coupling: a quantitative change - speed- 
up of He burning - for these stars leads to a qualitative 
change in the evolution - elimination of the blue loop. 
The elimination of the loop would have at least two ob- 
vious signatures, (i) An entire observed population of 
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stars, blue He burning giants, would not be accounted 
for. Detailed observations of blue loop populations exist 
(see, e.g., [56] and [57], particularly Fig. 3 therein), (ii) 
As stars go through a blue loop, they cross the Insta- 
bility Strip and become Cepheid variables. Without the 
blue loop [S3], one cannot account for the existence of 
Cepheid stars with the broad range of pulsation periods 
(corrseponding to ~ 8 — HA/ ). The initial crossing of 
this strip, as the star adjusts from its main sequence con- 
figuration to a He burning red giant state, proceeds too 
fast to give large enough numbers of these variable stars. 

Our investigations so far show that the resulting bound 
is somewhere between a rather conservative <?io < 0.8 and 
most aggressive g\$ < 0.5. The exact value depends on 
the detailed analysis of the observations and the physics 
of the simulation. While such a detailed study is well be- 
yond the scope of the present letter, below we summarize 
several relevant considerations. 

First, for our bound we require a complete disap- 
pearance of the blue loop, eliminating the entire ob- 
served population of stars. This is a conservative re- 
quirement. Given accurate counts, it may be possible 
to check whether the number of stars in the blue loop 
phase is reduced. For example, in Fig. [3] the middle panel 
shows that gio = 0.6 would reduce the time a 9.5Af 
star spends on the blue loop by a factor of two. (Notice, 
for comparison, that to get the same sensitivity for g\Q 
from solar-mass stars requires knowing the numbers of 
HB stars to a ~ 10% precision [40 .) 

Second, one can consider the effect of the axion on 
stars of different masses and find which stars have the 
most sensitivity. In our investigations, for example, we 
found that for 1O.5M stars the blue loop disappeared 
already for gio ~ 0.5. The observational signature in this 
case could be a gap in the observed periods of Cepheid 
stars, which vary as a function of stellar mass [BD] • Again, 
our bound is conservative with respect to this point. 

Third, the details of the blue loop depend on the treat- 
ment of the convection physics in the code [SS] [ST]. In 
our investigations with MESA we confirm that varying, 
e.g., the mixing length parameter shifts the exact value 
of g w at which the blue loop disappears. Understanding 
stellar convection is presently a focus of an active effort 
in the stellar astrophysics community. Since our code is 
being made public, we invite the members of this commu- 
nity to test the impact of various convection prescriptions 

- and other physical assumptions and numerical methods 

- on the axion bound. We hope, in time, this will result 
in a stronger bound on the axion. Tentatively, here we 
choose to state the conservative bound, gio < 0.8. 

CONCLUSIONS 

We have obtained a new astrophysical bound on the 
axion-photon interactions, by considering the evolution 



of stars ^7—12 times more massive than the Sun. The 
sensitivity of these stars to the axion-photon coupling 
compares favorably to the published bound gio < 1 from 
the solar mass stars [33] [55]. Sufficiently large axion- 
photon coupling is shown to eliminate the blue loop stage 
of the evolution, leaving one without an explanation for 
the existence of Cepheid stars in a broad range of pulsa- 
tion periods. This is the second time massive stars are 
used to constrain particle physics beyond the Standard 
Model and, as in the case of neutrino magnetic moment 
[43] . axion is also capable of qualitatively changing the 
stellar evolution. 
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FIG. 3: Impact of the A77 coupling on the evolution of a 9.5 M Q star through the He burning stage. Shown are the cases 
910 = {top), gio — 0.6 (middle), and gio = 0.8 (bottom). The horizontal axis shows the age of the star, in 10 6 years. Axion 
losses speed up He burning, resulting in the shortening and eventual elimination of the blue loop stage. 



